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We report the results on heavy quarkonium production in p + p and Au+Au collisions at midra- 
pidity via the dielectron decay channel at y^vJV = 200 GeV from the STAR experiment at RHIC. 
Results on the J/y pj spectra in p +p collisions are presented for 2 < pj < 10 GeV/c. The 
B-meson feed-down contribution to the inclusive J/y yield has been obtained using the J/y- 
hadron azimuthal angular correlation in p+p collisions, and is found to be 10-25% in the 
range 4 < pj < 12 GeV/c. The pj spectrum and nuclear modification factor for J/y with 
Pt < 10 GeV/c is reported, along with results from T(1S+2S+3S) production in Au+Au col- 
lisions. The nuclear modification factor for high-/?7- J/y (pj > 5 GeV/c) is found to be consis- 
tent with unity in peripheral collisions, while a significant suppression of low-pr J/y and T is 
observed in central collisions. The elliptic flow of J/y is reported for 20-60% central Au+Au 
collisions, and is found to be consistent with zero. 
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1. Introduction 

The suppression of J/y and other heavy quarkonium is expected due to the Debye color- 
screening of the potential between heavy quarks in a dense medium, and has been suggested 
as a signature of the formation of Quark Gluon Plasma (QGP) [1]. A suppression of J/y pro- 
duction in heavy ion collisions has been observed by the NA50 and NA60 experiments at the 
CERN-SPS [2] [3]. A similar amount of J/y suppression was observed at RHIC at midrapidity 
in y/s~^ = 200 GeV Au+Au collisions [4], despite the increased collision energy as compared to 
SPS. This can be understood by considering other effects due to the presence of a QGP which 
may contribute to the modification of heavy quarkonium production, such as statistical coales- 
cence of heavy quark-antiquark pairs or co-mover absorption [5, 6]. There are also contributions 
to the inclusive yields from feed-down effects [8] and the sequential melting of excited J/y and 
T states [7]. Cold Nuclear Matter (CNM) effects [9], such as the modification of nuclear PDFs 
(shadowing) [10], and final state nuclear absorption [11], need to be taken into account in order to 
fully quantify an anomalous suppression. This can be achieved by studying the production of var- 
ious quarkonium states in p +p, d+A and A+A collisions. Furthermore, p+p collisions can offer 
insight to the quarkonium production mechanism, as no model can yet fully explain the observed 
J/y yields in elementary collisions. The J/y elliptic flow (V2) can provide information about the 
contribution from recombination of charm quarks and the degree of thermalization of charm quarks 
in the medium. 

The results for the J/y pj spectrum in p + p and Au+Au collisions at ^fs^ = 200 GeV 
recorded by the STAR detector in 2009 and 2010 are presented in this paper. The centrality depen- 
dence of the J / y and T(1S+2S+3S) nuclear modification factor (R^a) are also shown. The elliptic 
flow of J/y is reported for 20-60% central Au+Au collisions at ^/s^ = 200 GeV. 



The reconstruction of J/ y and T are done via the dielectron decay channel, J / y, T — > e + + e , 
with a branching ratio (B) of 5.9% for J/y and 2.5%, 1.9%, and 2.2% for the IS, 2S, and 3S T 
states, respectively. The reconstruction methods used in this analysis are similar to those used in 
the 2005 and 2006 data [12, 13], and are described in [14, 15, 16]. The high-p^ J/y results shown 
in this paper are from 2009 p + p collisions with an integrated luminosity of l.&pb~ l , and 2010 
Au+Au collisions with an integrated luminosity of l.4nb~ l . The data were collected using the Bar- 
rel Electromagnetic Calorimeter (BEMC) to trigger events with a single tower (Atj x A0 = 0.05 x 
0.05) above the transverse energy threshold Ej > 2.6 GeV in p + p collisions and Ej > 4.3 GeV 
in Au+Au collisions. Over 300M minimum bias and 250M central-biased Au+Au events were also 
recorded, which allow for a precise low-py (pj < 5 GeV/c) J/y analysis. The Time Of Flight 
(TOF) detector, which covers full azimuth and \r\ \ < 0.9, was installed in 2009 with 72% com- 
pleted for recording p+p data, and was fully installed in 2010 for the Au+Au data recording. The 
Time Projection Chamber (TPC), which is responsible for momentum and dE/dx reconstruction, 
has been used along with the TOF and BEMC to obtain a J/y signal with increased signal-to- 
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background ratio and statistical precision. 
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Figure 1: The J/y pr spectra (left) and B-meson feed-down contribution (right) in p+p collisions. 



The left panel of Fig. 1 shows the corrected pj spectrum for J/y in p+p collisions at 
xfs,^ = 200 GeV (closed square). This is consistent with the previous STAR (downward trian- 
gle) [12] and PHENIX (open square) [17] measurements. A Blast-Wave model based on Tsallis 
statistics (TBW) [18, 19] has been fitted to all data points (solid line). Predictions from NRQCD us- 
ing LO color-singlet (CS) and color-octet (CO) [20] transitions are shown (dashed line), and agree 
with the data reasonably well but leave no room for feed-down effects from B, % c , and y/, which 
may account for up to 50% on the inclusive yield. The predictions from NNLO* CS [21] (grey 
band) underestimate the yield and pj shape at high-/^. The Color Evaporation Model (CEM) has 
been shown (dot-dashed line), and can explain the data reasonably well [22]. The right panel of 
Fig. 1 shows the pj dependence of the B -+ J/y feed-down contribution (closed circle). These 
results are extracted from J/i//-hadron azimuthal angular correlations using PYTHIA, and indicate 
a feed-down contribution of 10-25%, which is similar to measurements at higher collision ener- 
gies [23, 24, 25, 26, 27, 28]. 

The corrected pj spectrum for low-p^ J/y (solid black symbols) and high-pr (open black 
symbols) are shown in Fig. 2 (left panel), and are consistent with the previously published data 
(open red symbol) [4]. A TBW prediction using the freeze-out properties from light hadrons 
(dashed line) is compared to the data, and agrees with the data for pj > 2 GeV/c. The agree- 
ment is improved when fixing the radial flow (j8) to zero. This may suggest that the J/y has a 
smaller radial flow than light hadrons, or that regeneration from charm quarks may be a significant 
contribution at low-pr. 

The right panel of Fig. 2 shows the nuclear modification factor of J/y and T(1S+2S+3S) as a 
function of the number of participant nucleons (./Vpart) in Au+Au collisions at y/s^ = 200 GeV. 
The low-pr (2 < pj < 5 GeV/c) J/y (open star) and T (upward triangle) show a suppression that 
increases with increasing collision centrality, with Rj^a = 0-3 — 0.4 in central collisions. These 
results are consistent with previous measurements (open circle) [4]. The high-py {pj > 5 GeV/c) 
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Figure 2: The pr spectra for J/y (left) and nuclear modification factor for J / y and T(1S+2S+3S) (right) in 
Au+Au collisions. 



J/y is systematically above the low-px measurements for all centralities, and is consistent 
with unity in peripheral collisions (20-60%). This is consistent with previous high-p^ Cu+Cu mea- 
surements made by STAR (closed square) [12]. 
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Figure 3: The J/y elliptic flow, vi for 20-60% central collisions. 



Primordial J/ys which are produced in the initial hard scattering are expected to carry very 
little flow, however those which are subsequently created from the statistical coalescence of charm 
quarks may carry non-zero elliptic flow. The J/y elliptic flow is shown in Fig. 3 for mid-central 
Au+Au collisions at ^/s/^ = 200 GeV/c (closed circle), and is consistent with zero within errors. 
This is compared to the V2 from charged hadrons (open circle) and -mesons (upward triangle), 
which exhibit a non-zero elliptic flow. The small V2 of J/y in mid-central collisions may indicate 
that the complete thermalization of charm quarks is not achieved, or that recombination of charm 
quarks is a small contribution in this kinematic range. 
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